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ABSTRACT 

Active galaxies hosting two accreting and merging super-massive black holes (SMBHs) - dual Active 
Galactic Nuclei (AGN) - are predicted by many current and popular models of black hole-galaxy co- 
evolution. We present here the results of a program that has identified a set of probable dual ACN 
candidates based on near Infra-red (NIR) Laser Guide-Star Adaptive Optics (LGS AO) imaging with 
the Keck II telescope. These candidates are selected from a complete sample of radio-quiet Quasi- 
stellar Objects (QSOs) drawn from the Sloan Digital Sky Survey (SDSS), which show double-peaked 
narrow AGN emission lines. Of the twelve AGNs imaged, we find six with double galaxy structure, 
of which four are in galaxy mergers. We measure the ionization of the two velocity components in 
the narrow AGN lines to test the hypothesis that both velocity components come from an active 
nucleus. The combination of a well-defined parent sample and high-quality imaging allows us to place 
constraints on the fraction of SDSS QSOs that host dual accreting black holes separated on kiloparsec 
(kpc) scales: ^ 0.3 — 0.65%. We derive from this fraction the time spent in a QSO phase during a 
typical merger and find a value that is much lower than estimates that arise from QSO space densities 
and galaxy merger statistics. We discuss possible reasons for this difference. Finally, we compare the 
SMBH mass distributions of single and dual AGN and find little difference between the two within 
the limited statistics of our program, hinting that most SMBH growth happens in the later stages of 
a merger process. 



1. INTRODUCTION 

Active Galactic Nuclei (AGN) fueled by major galaxy 
mergers are predicted from models of black hole-galaxy 
co-evolution, in which black holes and their host galaxies 
grow together as an outcom e of hierarchical structure 
formation in the unive rse (e.g. lFerrarese k, Merrittll2Q00l : 
IGebhardt et"all [2000t) . In this scenario, powerful tidal 
forces funnel material into the centers of merging gas- 
rich galaxies, driving the growth of super-massive black 
holes (SMBHs) at near-Eddington rates. One prediction 
of this type of model is that prior to the final merger 
event, some ongoing gas- rich mergers should host dual 
AGNs separated by spatial scales comparable to the sizes 
of the galaxies. 

Two sets of factors determine the fraction of merg- 
ing galaxies that harbor dual AGN. One is the causal 
link between galaxy mergers and high accretion rate 
phases of SMBH growth. Studies show that most low- 
luminosity AGN s (i.e., Seyferts) are n ot undergoing sig- 
nificant mergers (|Cisternas et al.ir2010( ) , though they may 
be frequently associa ted with galaxy- scale tidal interac- 
tions revealed in HI ()Kuo et al.ll2008[ ). However, power- 
ful AGNs (i.e. QSOs) are believe d to be strongly cor- 
related with major n i erger events ([Hopkins et al.ll2006l : 
iHopkins fc HernquistI 120091 ). Such studies have been 
borne out by associations between statistically deter- 
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mined quantities such as the quasar luminosity function, 
galaxy and merger rn ass functions, and quasar clustering 
(jHopkins et"all[2008l ). 

The second key factor is the frequency at which both 
black holes in a merger accrete gas and shine together as 
AGN, while still separated on scales of a few kpc. This 
fraction is not well constrained. Most computational 
models do not address pre-coalescence fueling very accu- 
rately, since it depends on several parameters that lie be- 
low the resolution of typical simulations. If a substantial 
fraction of QSOs show dual active nuclei and signatures 
of strong activity in the early or intermediate stages of a 
galaxy merger, a significant portion of the growth of the 
SMBH will occur before final coalescence. On the other 
hand, if a great majority of QSOs are seen as single AGN, 
the SMBHs in these systems have either in-spiralled to 
parsec separations or have already coalesced before sub- 
stantial activity begins. In such a case, most black hole 
growth is expected to occur in the late part of a merger. 
Knowledge of the phase during which most of the accre- 
tion and QSO activity takes place has important conse- 
quences for the form and scatter of SMBH scaling laws, 
because black hole growth in a gas-rich galaxy merger 
runs in parallel with another important avenue of gas 
depletion: star-formation. If a large fraction of SMBH 
growth occurs early in the merger, the eventual relation- 
ship between the mass of the black hole and the stellar 
content of the merger remnant will be influenced. 

The incidence of dual AGNs in an ongoing merger 
therefore provides insight into the physical processes tak- 
ing place when the two merging SMBHs are separated 
on galaxy scales or less. In the past few years, studies 
have searched large extra-galactic spectroscopic datasets 
to arrive at statistically useful samples of double QSO 
or AGN candidates ()Hennawi et al.. .2006; .Mvcrs et al.l 
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120071 120081: iHennawi et"all I2010D with separations on 
the sky up to 650 kpc. More recent studies identi- 
fied pairs of Q SOs separated on scales of t ens of kpc 
(iForeman et a l. 2009) to a few kpc (Junkk arinen et al. 
200lVConierford ct al. 20"ol !Liuet al.il2009i : ISmith et al. 
2010,: Wang ct al. 2009). From such studies, it is clear 
that the fraction of candidate dual AGNs among QSOs 
is small, even after accounting for the short duty cy- 
cle of phases with high accretion rates and the relative 
time-scales during which both nuclei are visible as AGNs 
(jSmith et al. 2010). Note that the identification of dou- 
ble AGN separated by a few kpc from spectroscopic sur- 
veys gives an upper limit to the dual AGN fraction, since 
such studies select for kinematic substructure in narrow 
emission lines which may also arise from AGN jets or 
outflows. 

In this paper we approach the issue by looking for the 
"smoking gun" : individual galaxies which contain dual 
AGNs and whose morphologies are typical of merging 
systems. To place constraints on the fraction of QSOs 
which show double structure on scales of 10 kpc or less, 
we have undertaken a program of near-IR (NIR) imag- 
ing of candidate radio-quiet double QSOs using the Keck 
Laser-Guide Star Adaptive Optics (LGS AO) system 
(jWizinowich et al][2006i ). In this Paper, from our 12 tar- 
gets we report six dual AGN candidates, four in close 
mergers and two in double systems, with separations 
ranging from 3 to 12 kpc. We measure luminosities of the 
AGN and host galaxies, the separations of the merging 
components, as well as black hole masses of the brighter 
QSO in each pair. Since we start with a complete par- 
ent sample of QSOs, we are able to refine the fraction 
which show double structure on scales larger than 0.45- 
1.6 kpc (over our redshift range), and we discuss the 
frequency of such closely separated merging AGN in the 
SDSS QSO population at large. We com pare our results 
with a recent study by iFu et al.l (|2010D which images a 
larger sample of Type I and Type II AGN to a shallower 
depth. Many of the doubl e systems presen ted in this pa- 
per can also be found in iFu et al.1 (j2010f) . We adopt a 
Concordance Cosmology with Hq ~ 70 km/s/Mpc and 
flA = 0.7. 

2. KECK LASER GUIDE STAR AO IMAGING PROGRAM 

For a well-defined sample of pre-coalescen ce dual AGN 
candidates, we drew upon the catalog of Smit h et al.l 
(|201Ci) . which selected SDSS DR7 spectroscopic QSOs 
that show a double-peaked [O III]A5007 line. Many 
of them also show similar line-splitting in other bright 
emission lines, such as [O II] A3 72 7 and II/3. The 
full catalog contains 86 Type I AGN. We excluded 
radio-detected QSOs, which are frequently associated 
with disturbed emission liiie regions due to jet interac- 
tions (jWilson fc WilUsllT9MIWhittiaiT99llR,osario et al.l 
I20T0I) . We also excluded Type II AGN in the SDSS 
QSO database which are frequently failures in the QSO 
spectral identification algorithm and hence have rather 
unconstrained selection criteria (jSmith et al.1 '2010). To 
our working sample, we applied LGS AO observability 
requirementf0 which yielded 37 possible LGS AO tar- 
gets. These were then ranked according to the brightness 

^ An r < 17.5 tip-tilt (TT) star within an arcminutc of the 
galaxy for a good tip-tilt correction. 



and nearness of their tip-tilt (TT) stars (which deter- 
mines the performance of the AO system), the strength of 
galaxy spectral features in their SDSS spectra (which de- 
termines the contrast between the AGN point source and 
the extended galaxy light), their observed SDSS magni- 
tude (which determines the S /N of the final images) and 
the quality of the [O III]A 5007 line spli t , follo wing criteria 
discussed and applied in ISmith et al.l ()2010[ ). 

Over a full night on Sep 20 2009 and two half-nights on 
Jan 4 and Apr 20 2010, a total of 12 QSOs were imaged. 
The NIRC2 near-infrared camera (PI: Keith Matthews) 
on the Keck II telescope was used in the Wide-Camera 
mode (FOV: 40", pixel scale: 0.04") with a broadband H 
filter. After imaging each object, a Point-Spread Func- 
tion (PSF) reference star pair was imaged with similar 
TT star properties as those of the science target. This 
ensured that the PSF reference suffered consistent aniso- 
planatic effects as the target QSO. Where a good PSF 
reference pair was not available, a series of short expo- 
sures of the TT star was taken instead. 

Exposure times ranged from 20 — 60 min depending on 
the brightness of the target. Large numbers of dithered 
sub-exposures were taken to prevent saturation of the sky 
background, and to overcome flat-field variations across 
the detector. A similar process was used for the PSF ref- 
erence pair images, though the sub-exposure integration 
times were much shorter to prevent detector saturation. 
Sub-exposures were combined in post-processing using a 
version of the NIRC2 IDL pipeline, modified to robustly 
account for the very narrow core of the PSF. This is crit- 
ical when modeling the QSO host galaxy to estimate the 
brightness of the central point source. 

3. PHOTOMETRY 

Photometry was performed directly from the NIRC2 
images. The NIRC2 detector in our default operating 
mode (Multiple-Correlated Double Sampling or MODS) 
is known to show a considerable excess of pixels with low 
counts, believed to be due to the rounding of low valued 
pixels at readout0. The excess tends to skew estimates of 
the local background from modal statistics. We instead 
fit a gaussian to the local pixel distribution around the 
object, excluding the pixel values around zero and use 
the mean of the gaussian as the local background. A 
typical absolute difference in the local sky estimate from 
modal statistics vs. the gaussian fit are around 20% but 
this can be more pronounced in images where the local 
background is higher. The effect of the low pixel excess 
on photometry after sky subtraction is negligible. 

All photometry was referred to the 2MASS H-band 
magnitude system. For objects which had a nearby 
bright star in the field, the star was used as a photo- 
metric calibration standard to the 2MASS system. For 
the rest, the total photometry of the QSOs themselves, 
measured in a large aperture (5" in diameter), was com- 
pared to measurements from the 2MASS point source 
catalog. The latter calibrations have a larger error, since 
most of the QSOs are barely detected in 2MASS. 

4. DOUBLE STRUCTURE IN THE H-BAND IMAGES 

The final H-band images of the 12 QSOs are presented 
in Fig. [TJ Our images resolve structure on scales of a kpc 

•5 See the NIRC2 News page 
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TABLE 1 

Sample of QSOs with NIRC-2 LGS AO imaging 



bUbb IJJ 


z 


Uoubic : 


Sep. ( ) 


bep. (kpcj 


m]Lf (QbUJ 


mpq- (host) 


m]Lf (companion) 


J153231. 80+420342. 7 


0.210 


n 


N/A 


N/A 


16.81 


16.54 


— 


J091649. 41+000031. 5 


0.222 


n 


N/A 


N/A 








J161027.41+130806.8 


0.229 


y 


2.35 


8.55 


16.87 


15.77 


19.98 


J081542.53+063522.9 


0.244 


n 


N/A 


N/A 


16.46 


16.01 




J095207.62+255257.2 


0.339 


y 


1.00 


4.82 


17.96 


16.36 


17.43 


J130724.08+460400.9 


0.353 


y 


2.37 


11.64 


20.47 


17.62 


17.86 


J020011. 52-093126.1 


0.361 


y 


1.17 


5.84 


16.31 


16.93 


19.55 


J140923. 51-012430.5 


0.405 


n 


N/A 


N/A 


16.65 


19.13 




J124859. 72-025730.7 


0.487 


y 


0.53 


3.15 


16.58 


17.22 


16.82 


J154107.81+203608.8 


0.508 


y 


2.00 


12.22 


17.88 


18.68 


19.87 


J121911. 16+042905. 9 


0.555 


n 


N/A 


N/A 


16.01 


20.10 




J072554.42+374436.9 


0.634 


n 


N/A 


N/A 


16.59 


18.74 





or less in all these galaxies. A large fraction of our tar- 
gets appear to have companions or show signs of being 
in a galaxy merger, highlighting the value of kinematic 
methods to find AGN in mergers. Of the twelve, three 
are in pairs with fairly equal H-band luminosity ratio: 
J095207+255257, J130724+460400 and J124859-025730. 
The obvious signatures of an interaction in these sys- 
tems strongly suggest that they are undergoing a merger 
event. J020011-093126 has a minor companion, but a 
low surface brightness bridge connects the two galax- 
ies. We will refer to these systems as 'mergers' in the 
rest of the paper. Two more galaxies appear to have 
significantly fainter companions: J161027+130806 and 
J154107+203608. These may be minor mergers or satel- 
lite systems, or the fainter object may be entirely unre- 
lated to the AGN host. We call these 'minor doubles'. 

For the systems with double structure, the separation 
between the centers of the two galaxy components are 
listed in Table 1, both in arcseconds and kpc. The 
four clear mergers have separations ranging from 12 kpc 
down to 3 kpc, and three of the four are separated by 
less than 6 kpc. Note that the angular diameter of the 
SDSS spectroscopic fiber aperture (3") corresponds to 
distances of 5 — 15 kpc over the redshift range of our 
galaxies. For some of the pairs with large separations 
(> 10 kpc), it is not clear that sufficient nuclear light 
from the neighboring galaxy will have entered the SDSS 
fiber to produce double-peaked [O III] lines, in which case 
the line splitting may be produced by processes related 
to just one of the nuclei. This is particularly relevant for 
J072554+374436, which has a compact point-like neigh- 
bor that is 2.9" away (not shown in Fig. [ij. This sep- 
aration is about the diameter of the SDSS aperture. In 
this case, it is highly unlikely that the second component 
to the [O III]A5007 line came from the nearby object, 
which, based on its compact appearance, may be a faint 
foreground star. In general, the chance of random align- 
ment between one of the QSOs and a faint un-related 
galaxy in the foreground or background is quite small (a 
few percent), based on source counts of galaxi es in the 
H-band from NIR surveys (e.g. lFrith et alfcoOGl) . There- 
fore, most of the double galaxies in this NIRC2 imaging 
sample are expected to be real associations. 

We modeled the H-band images of the AGNs usi ng the 
galax y structure fitting code GALFIT v3.0 ( Peng et al.l 
120021 ). The brighter galaxy was modeled as a combi- 



nation of a point source and an extended stellar light 
distribution having an elliptical Sersic profile with a vari- 
able Sersic index, half-light radius, ellipticity and orien- 
tation. In the cases with double structure, we modeled 
the fainter galaxy with a single elliptical Sersic profile, 
making the assumption that the point-like QSO is to be 
found only in the brighter object of the pair (as discussed 
in §5.2). Each fit was examined by eye, and obvious fail- 
ures were re-run with a different set of initial parameters 
until a good fit was achieved. The primary aim of these 
fits was to estimate the total H-band luminosity of the 
nuclear point source and the host galaxy, as well as the 
companion, if applicable. These estimates are listed in 
the last three columns of Table 1, where mjj is the ob- 
served mag nitude in the H-band in the 2MASS magni- 
tude system. In a forthcoming paper, we will present a 
more detailed and accurate treatment of the structural 
analysis of our sample, including comparisons to NIRC2 
images of inactive galaxi es from the Cente r for Adaptive 
Optics Treasury Survey (|Koo et al. II2007I ). 

While these targets are all selected to be optical Type 
I AGN, it is clear that not all are dominated by a central 
point source. Extended host light is seen in most of these 
galaxies. This is partly due to the fact that we chose to 
image, when possible, objects that appeared to have a 
non-negligible fraction of their light coming from the host 
in the SDSS spectra. The AGN point source contributes 
from 25% to 100% of the total light in most of these sys- 
tems, though in two cases, the estimated AGN fraction is 
much lower. The first, J130724.08+460400.9, is a pair of 
clearly interacting widely separated bright galaxies and 
the point source contributes less than 10% of the total 
light of the larger galaxy. The image of the second case, 
J091649. 41+000031. 5, shows an elliptical galaxy with a 
fairly smooth core and the GALFIT modeling suggests a 
negligible amount of light from a point source. However, 
our PSF monitoring images, taken just after the images 
of the galaxy, show that the performance of the AO sys- 
tem suffered significant variation during this phase of the 
night. A broad PSF, probably due to an unsettled low 
bandwidth wavefront sensor in the AO system, would 
smooth out the core of the image and lower the contrast 
of a point source. Therefore, we are unable to accurately 
estimate the point source fraction in this galaxy, but can 
conclude that it does not dominate the H-band light. 
This example serves to demonstrate the importance of 
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J0952+2552 




J1248-0257 




Fig. 1. — Keck laser guide star AO H-band images of double-peaked narrow line SDSS QSOs, taken using the NIRC2 camera. The scale 
bar corresponds to 4 kpc at the redshift of the system. Six of the targets show signs of merger activity or close companions. 



concurrent PSF monitoring in AO studies of the struc- 
ture of distant galaxies. 

5. DO BOTH GALAXIES HOST AGNS? 

We have shown that double galaxy structure is seen 
in a large fraction of QSOs which show double peaked 
narrow-lines. However, the question arises whether the 
fainter of the two galaxies also hosts an active nucleus, 
i.e, are these systems true double AGN? While sensitive 
to mergers, our NIR images do not provide color informa- 
tion about the galaxies or their point sources. Prom the 
single-band images alone, we are unable to test whether 
an AGN exists in each nucleus in the cases of the galax- 
ies with double structure. However, we can carry out a 
simple test which compares the ionization level of both 
[O III] velocity components seen in the SDSS spectra. 
This allows us to determine if one of the peaks comes 
from, say, a low metallicity star-forming satellite galaxy. 

We concentrate on the ionization-sensitive 
[O III]A5007/H/3 ratio which involves two lines that 
are close in wavelength and therefore are not strongly 
influenced by differential extinction or by the varying 
resolution of the SDSS spectrograph across its spectral 
range. This ratio was measured by first modeling the 
[O III]A5007 line from each object as the sum of two 



independent gaussian profiles, which were then scaled 
in flux to match the narrow H/3 line. In all cases, 
the broad-line H/3 and the local continuum level were 
concurrently fit by eye using a high order polynomial 
with a variable index. Uncertainties in the fitting of 
the continuum were incorporated into the error in the 
normalization of the gaussian components. 

In Fig.H we show [O III]A5007/H/3 for each of the two 
[O III] lines in the twelve QSOs in our sample. From this 
figure, it is clear that most of our measured line ratios are 
above the highest possible star-forming value for normal 
galaxies, so these components are most likely AGN. Al- 
most all lie to the right of the middle star-forming value, 
implying that AGN ionization dominates the narrow H/3 
line in most of these objects. 

The one case with low [O III]A5007/H/3 ratios, 
J020011. 52-093126.1, is surprising, since the image 
clearly shows a bright point source and the spectrum 
is dominated by strong QSO spectral features. There is 
no doubt that a bright AGN is present in this system. 
Interestingly, this is a double system showing some sig- 
natures of a merger. Both velocity peaks show a low 
ratio, indicating that the H/3 line is enhanced relative 
to [O III]A5007, with respect to the typical AGN, for 
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Fig. 2. — A graphical comparison of [O III]A5007/H^ ratio of 
the two velocity components, measured from SDSS spectra, of the 
QSOs with NIRC2 imaging. Imaged objects with no double struc- 
ture are shown with black points, while red is used for objects with 
large luminosity ratios between the galaxies in the image (minor 
doubles) and green is used for objects with roughly equal luminos- 
ity ratios (mergers). The best estimates for the ratio are shown 
as filled points, while the 2a lower limits to the ratio are plotted 
with open points; dotted lines join the measured values and limits, 
to guide the eye. Two objects have undetected narrow H/3 and 
therefore only have lower limits on the ratio. The vertical dashed 
lines are three representative values of [O III] A5007/H/3 along the 
curve separating star-forming galaxies from composite/AGN sys- 
tems, as defined by Kauffmann et al. (2003). The highest value 
(~ 8) is for the lowest mass star-forming galaxies generally found 
in the SDSS, the middle value of ~ 5 is for typical low-mass star- 
forming galaxies and the lowest value corresponds to typical high 
mass galaxies. The majority of the velocity components have ratios 
consistent with AGN ionized gas. 

both components in this system. The low level of ioniza- 
tion in both peaks suggests that both components come 
from a common ionized region, supporting an outflow as 
a likely explanation for the line-splitting, rather than a 
merger. This implies that some QSOs that are in merg- 
ers may have a double [O III] profile that is unrelated to 
the merger activity. 

Five of the AGNs with double structure in our im- 
ages lie to the right of the vertical line for which 
[O III]A5007/H/3 = 8, whereas only two of the non-double 
AGNs are in this region. This may point to a subtle trend 
that AGN in close double systems tend to have higher 
overall levels of nuclear ionization, possibly due to dif- 
ferent nuclear gas environments mediated by the inter- 
action. However, a significantly larger sample of isolated 
and merging AGN is needed to test this difference in a 
statistically meaningful way. 

6. DISCUSSION 

6.1. AGN luminosities 

In Fig. 131 we compare the total 5100 A monochro- 
matic luminosity distributions of the 12 AGNs in our 
sample (filled histogram) with that of the parent sample 
of double-peaked Type I AGN from the SDSS (solid open 
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Fig. 3. — Distributions of monochromatic luminosities at 5100 A. 
The open solid histogram is the distribution for the parent sam- 
ple of SDSS Type I AGN with double-peaked narrow lines, which 
may be compared to the filled histogram, which is the distribution 
for the 12 AGNs with H-band imaging presented in this paper. 
The dashed histogram is the distribution for 28 PG-QSOs from 
the QUEST survey. Our sample is equally divided among high 
luminosity Seyferts and low luminosity QSOs. 

histogram) . The traditional divide between Seyferts and 
QSOs (log J^Lsioo — 44.37) is indicated by a vertical 
dashed line. Our sample has a range of about an or- 
der of magnitude in luminosity and spans the divide 
between luminous Seyfert Is and QSOs, with a median 
value of 44.32, slightly below the median of the parent 
sample. We also include in the figure, as a dashed his- 
togram, the distribution of luminosities for a sample of 
local PG-QSO s from the Quasar/U LIRG Evolution Sur- 
vey f QUEST. IVeilleux et 3.1.1 1200^). Our objects overlap 
with the faint end of the QUEST distribution, but we 
do not have any AGNs that sample the high luminos- 
ity end among local QSOs, which can be more than an 
order of magnitude brighter than our brightest objects. 
One must keep in mind that the following discussion is 
valid only for AGNs that lie at the transition between 
normal Seyferts and QSOs, and our results may not be 
applicable to the most luminous QSOs. 

6.2. Black Hole Masses and Growth 

With a set of genuine Type I AGN in pre-coalescence 
mergers, we can estimate SMBH masses relative to the 
masses of the merging hosts and may be able to test, 
with a sufficiently large sample, whether the black hole 
growth leads or lags the growth of stellar mass in the 
hosts. While our current sample of double AGN is too 
small to make many statistically significant statements, 
we nevertheless explore this possibility in order to under- 
stand the limitations and complexities of such an analy- 
sis. 

We start by making the weak assumption that each 
QSO system contain only one Broad-Line AGN. The sec- 
ond active nucleus, if present, is a Type II AGN, i.e, ob- 
scured along our line of sight so that only the narrow-line 
emission is visible. Surveys of local AGN find that Type 
lis outnumber Type Is by a factor of 3, possibly deter- 
mined by the mean covering factor of the putative dust 
torus believed to exist in most AGN (Risaliti ct al...l999.) . 
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Fig. 4. — Black Hole masses of the double-peaked [O III] QSO sample. The left panel (a) compares the distribution of Mbh of the parent 
SDSS double- [O III] sample (open histogram, scaled down in number) with the masses of all 12 QSOs in our sample (solid histogram), 
uncorrected for host galaxy contamination. The galaxies in our sample typically host more massive black holes than the average SDSS 
double-[0 III] QSO. The right panel (b) compares Mbh between AGN in double systems (solid line) and isolated systems (dashed line), 
corrected for the light of the host galaxy. While a difference in the distribution is seen, it is not significant, given the small number of 
objects in each sub-sample. 



However, there is some evidence that thi s difference ma y 
be smaller for distant or luminous AGN ()Hasingeiil2008[ ) . 
Since our dual AGN are typically at the boundary be- 
tween Seyferts and QSOs in terms of their nuclear bolo- 
metric luminosity, we will assume that Type II nuclei 
are more common and that all of the dual AGN in our 

sample are pairs of Type I and Type II nucl ei. 

We follow the prescription of lShields et alj ()2003l Eqn. 
2) and estimate the masses of the SMBHs as follows: 

Mbh = 10'' ''^^'3ooo-^44' solar masses (1) 

where -uaooo is the FWHM of broad in units of 3000 
km s~-^ measured from the SDSS spectra and L44 is the 
continuum luminosity of the QSO per decade in wave- 
length at a rest wavelength of 5100 A, in units of 10_f^_erg 
s~^. W e follow the technique outlined in Salviandc r et alJ 
(|2007| ) for the measurement of the H/3 line shape from 
SDSS spectra. After subtracting the narrow H/3 emis- 
sion line using a 1:10 scaled version of the [O III]A5007 
line, the remaining broad H/? line is modeled as a Gauss- 
Hermite function, from which W3000 is derived. 

The continuum luminosity L44 is also measured di- 
rectly from the SDSS spectra of the QSOs. The light 
in the large aperture SDSS spectrum is a combination of 
emission from the QSO and from the host galaxies in the 
system. Therefore, we need to apply a correction factor 
to which accounts for host galaxy light within the 
SDSS spectral aperture. The correction term was esti- 
mated in the following way. First we calculated the flux 
ratio between the observed H-band and the rest-frame 
5100 A (essentially a color: Hots ~ 5100) from the Spec- 
tral Energy Distributions (SEDs) of an average SDSS 
QSO and a set of normal galaxies, taken from the pub- 



licly available SWIRE template library p'olletta et alJ 
I2007D . We calculate this color at the redshifts of each of 
the QSOs in our sample. We found that the difference in 
Hobs — 5100 for early-type and late-type galaxies is small 
at all our redshifts (< 0.5 mag), so we adopted a typical 
value. 

We then smoothed our NIRC2 images to approximate 
ground-based seeing and, guided by our GALFIT fits to 
our NIRC2 images, measured the flux in the point source 
relative to the total light within a 3" circular aperture 
centered on the galaxy (or the brighter galaxy, in case of 
doubles). This measurement is equivalent to determin- 
ing the point source fraction in an SDSS-like aperture, 
but in the observed H-band. Applying our estimates of 
Robs — 5100, we scaled the H-band point source fraction 
to get the point source fraction at rest-frame 5100 A. 
Since QSO SEDs are much bluer than galaxy SEDs, the 
point-source fraction increases from a modest fraction 
of the hght in the H-band (greater than 10%) to 40% 
to 100% of the light at 5100 A. Finally, we applied the 
point-source fraction at 5100 A to our estimates of ^44 
measured directly from the SDSS spectra, to get revised 
values of the QSO continuum luminosity. It transpired 
that, since Mbh depends only on the square- root of L44, 
the effect of this correction is only modest (less than 0.1 
dex reduction in the median M^^f of our sample). 

In FigH^, the distributions of the uncorrected black 
hole masses (i.e, before applying the point-source fraction 
correction to L44) of the QSOs in our sample (shaded his- 
togram) are compared t o the scaled d istrib ution of the 
parent SDSS sample of ISmith et al.l ()2010 f) (open his- 
togram). We have excluded J091649.41-h000031.5, as the 
point source fraction is very uncertain for this particular 



Double QSOs from SDSS 



7 



galaxy (see §3). The median uncorrected BH mass of 
our AGN (log Mbh = 8.43) is significantly larger than 
that of the parent sample (log Mbh = 7.94). This differ- 
ence is probably an effect of sample selection. Since we 
specifically targeted QSOs with some indication of host 
galaxy light in their SDSS spectra, we are likely to select 
systems with massive host galaxies and, by extension, 
massive black holes. Note that this same bias may come 
into play in studies that estimate QSO host g alaxy ve- 
locity dispersions from SDSS spectra, such as in lFu et aD 
(jlOlO). Therefore, care must be taken in the interpreta- 
tion of SMBH scaling relationships derived with these 
methods. 

In Fig|4l3, a comparison is made between the cor- 
rected Mbh distribution (i.e, after accounting for the 
host galaxy) of the six QSOs which are in mergers or have 
companions against the corrected distribution of the five 
isolated QSOs (again excluding J091649.41+000031.5). 
The median value of log M bh for the double systems is 
slightly lower than that of the isolated systems, by 0.14 
dex, but the difference is not significant ( as determined 
by a two-sided Kolmogorov-Smirnov test). There is no 
statistical difference in this result if we restrict ourselves 
to only the four objects with major mergers. Therefore, 
we conclude that we do not detect any major differences 
in the Mbh distribution between pre-merger and normal 
QSOs. If much SMBH growth occurs in QSOs, a sub- 
stantial portion may not occur in the pre-merger phase. 
Major caveats do remain, however. If, for example, a 
large fraction of such systems harbor dual Type I AGN, 
in contrast to our in-going assumption in this analysis, 
then the Mbh we calculated above are too high and will 
have to be adjusted. A larger sample of QSOs in mergers 
will help test this in the future. 

6.3. The Incidence of Double AGN 

Taken at face value, the presence of 6 double or merg- 
ing galaxies out of our sample of 12 QSOs implies that 
between 30% - 70% of Type I radio-undetected double- 
[O III] emitters are in galaxy pairs at 0.2 < z < .6, using 
binomial statistics of small numbers following iGehrelsl 
([T986,). Of these, we find 4 (16% - 56%) in major merg- 
ers. These fractions compare favor ably with thos e de- 
rived by a similar LGS AO program (iFu et al.ll2010l) over 
the same redshift range. Given the small number of ob- 
jects in our sample, we are unable to make any substan- 
tive statements about the evolution of this fraction, as 
iFu et al.l (|2010l l have done, though they use a broader 
parent sample, including radio-detected and radio-loud 
AGN. We caution that the visibility of faint companions 
to the QSO hosts is strongly affected by surface bright- 
ness dimming towards higher redshifts and therefore a 
full treatment of the evolution of this fraction needs to 
take these factors into account. 

We have verified, using H-band number counts from 
the literature, that chance alignments with background 
or foreground galaxies are not a likely explanation for 
the companions that we find around our QSOs. How- 
ever, massive galaxies, such as those that host QSOs, 
typically reside in dense galactic environments. There- 
fore, the possibility exists that some or all of the minor 
doubles may be due to nearby satellite galaxies rather 
than ongoing mergers. 

To estimate the fraction of normal QSOs with faint 



neighbors, we examined a set of HST/NICMOS images 
of 28 QSO hosts published as part of the QUEST survey. 
The median redshift of the QUEST QSOs is lower than 
ours (about 0.14, compared to 0.35 for our NIRC2 sam- 
ple). Despite this difference, we choose to use this com- 
parative sample since it is derived from a homogenous 
high resolution H-band imaging dataset of normal QSOs 
with high quality and depth. Since the QUEST QSOs 
are not pre-selected to be double-peaked AGN and since 
the vast majority of QSOs do not show double-peaked 
narrow lines, we expect that the QUEST QSOs are pre- 
domina ntly normal QSOs that would not have made it 
into the lSmith et al.l ()2010D catalog had they been at sim- 
ilar redshifts. Therefore, the neighbor statistics of these 
QSOs should be generally app licable to the SDSS QSO 
parent sample from which the iSmith et al.l ()2010l ) sam- 
ple was selected. We also assume that there has not been 
major evolution in the population of QSO host satellites 
between z ~ 0.1 and z ^ 0.4. 

From an examination of the QUEST images, we de- 
termine that approximately 35% of 22 radio-quiet QSOs 
have a fainter companion within 15 kpc, but only two 
(~ 9%) show obvious spatially-resolved signs of merger 
activity. To compare, our double-peaked QSO sample 
has roughly the same fraction of minor doubles, but a 
larger fraction of mergers. Under the assumption that 
the statistics for the QUEST sample apply to the SDSS 
parent population of QSOs, we conclude that some or all 
of the minor doubles in our sample may arise from satel- 
lite galaxies around QSO hosts. If these satellites also 
host AGN, then the double-peaked narrow lines may in- 
deed come from two separate active nuclei. However, it 
is also possible that the double-peaked structure in the 
minor doubles could come solely from the active nucleus 
containing the QSO, as in the 50% of the NIRC2 sample 
that do not show double structure. At present, we are 
unable to differentiate between these two possible expla- 
nations for the minor doubles. 

On the other hand, the rate of major mergers in our 
sample is considerably higher than those seen among 
QUEST QSOs, which implies that these mergers are 
likely to be the source of the double-peaked lines and are 
probably double AGN. In addition, the two cases of on- 
going massive mergers among the QUEST sample appear 
qualitatively different from the mergers in our sample, in 
that they seem to have double nuclei that are at closer 
separations and are probably in a more advanced merger 
stage. 

6.4. Fraction of QSOs with Dual SMBHs 

Assuming that each of the two [O III]A5007 peaks cor- 
responds to a galaxy in the double systems, we are able 
to calculate the fraction of QSOs that host dual active 
SMBH s on scales of few to several kpc. ISmith et al.l 
()2010[) estimate that 0.9% of radio-undetected QSOs 
show double-peaked narrow AGN lines. Applying the 
fraction that show resolved double structure in our 
NIRC2 images, about 1 of every 200 QSOs (0.5%) con- 
tain double accreting SMBHs. If, instead, we assume 
that the minor doubles are not double AGN, this ratio 
drops to 0.32%, or 1 in 300 QSOs, though this may be as 
high as 0.5% at the la upper uncertainty in our estimate. 
We denote this fraction by fduai- Note that we inde- 
pendently find the same fraction of mergers as iFu et all 
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(|2OT0l) among double-peaked AGN. Since lFu et all (|20lol ) 
present a larger sample of imaged objects, we expect that 
the total uncertainty in fduai should be smaller than our 
estimate above, but we will adopt the more conservative 
uncertainty from our limited sample. 

We can compare fduai to a rough estimate of the num- 
ber of dual AGN expected among QSOs. Given the 
assumption that all QSOs are associated with a major 
merger event, let F be the 'duty cycle' of QSO activity 
during a typical merger, i.e, it is the fraction of the time 
during a typical merger when either one of the SMBHs 
is active. Let A be the fraction of the time in which a 
typical merger is in its early stages, with nuclei separated 
by kpc distances or more. Then, the fraction of QSOs 
which have double nuclei that we may observe with our 
NIRC2 program is: 

fduai = TA (2) 

where we have assumed that the activity of the two 
SMBHs in a merger are independent of each other, and 
that the activity also does not depend on the stage of 
a merger. This last assumption is unlikely to be true - 
simulations suggest that most of the activity happens at 
the latter stages of a merger. We include it here for sim- 
plicity and use our final calculations to test its validity. 

Morphological studies of mergers at various stages 
of development, as w ell as numerical simulations (e.g., 
iCallegari et al.|[2009l) . imply that most merging galaxies 
spend a substantial fraction of the time of the merger in 
a pre-coalescence phase with widely separated nuclei. In 
other words, A should be a fairly large fraction. Given 
fduai ~ 0.003 — 0.0065, Eq. 2 places F in a similar range. 
How does this compare to simple estimates of the QSO 
merger 'duty cycle' from population statistics? 

If we assume the QSOs are found only among rela- 
tively massive galaxies (logM* > 10.5 Mq), we can place 
some constraints on F by comparing the merger rates of 
such galaxies with the frequency of QSOs among them. 
Based on the analysis of morphologies and pair statis- 
tics from field galaxy surveys, roughly 4% of massive 
M*) galaxies host on-goin g galaxy mergers dB undv et al.l 
12009; Hopkins ct alj|20100. This fraction has some uncer- 
tainty, but is unlikely to be greater than 10%. From the 
ratio of the spa ce densities of 10^^ M0 galaxies lO"'^-^ 
Mpc-3 dex-^ IbH^v et al.l 120061 ) and QSOs in the lu- 
minosity range of our sample (—23 < My < —21, 10^^ 
Mpc-3 mag~ MCroom et al.ll200l . at z - 0.5, QSOs are 
found in only about 0.3% of massive galaxies. This frac- 
tion is roughly consistent wit h the more detailed anal- 
ysis by iHamiiton et"al] (f200l at z = 0.26. The QSO 
'duty cycle' F is given by the ratio of the fraction of 
QSOs among galaxies to the fraction of mergers among 
galaxies, as long as all QSOs are associated with merg- 
ers. Therefore, F ~ 0.003/0.04 ^ 0.08. This estimate is 
significantly different than the value we derived earlier 
from dual AGN statistics. Put simply, the fraction of 
time a merging pair of galaxies spends in a QSO phase is 
more than an order of magnitude higher than an estimate 
based on our measurement of fduai ■ 

There are several possible explanations for this discrep- 
ancy and we briefly highlight a few here. The simplest is 
that our estimate of the value of fduai, with its consid- 
erable uncertainty, is not representative of the true pop- 



ulation of AGN in this luminosity range. How ever, our 
fractions match well with a larger sample from lFu et al.l 
([2010') and are unlikely to be incorrect by an order of 
magnitude (i.e, we measure only 30% in mergers, when 
the true fraction is almost 100%). 

An explanation that is outlined in iSmith et al.l (j2010l ) 
suggests that searches for double-peaked emitters are 
highly incomplete and miss as many as 90% of true merg- 
ers, due to line of line-of-sight effects. This would in- 
crease fduai and bring into closer agreement the two es- 
timates of F. Alternatively, a majority of SDSS Type I 
AGN may not be associated with massive mergers, but 
driven by processes like minor mergers or secular infall, 
as recently proposed for X-ray selected Seyfert galaxies 
with s imilar AGN lum inosities at intermediate redshifts 
(Ciste rnas et al. j[ 2010D. High reso lution studies of QSO 
host structure (jGuvon et al.l '2006) find that significant 
fractions (~ 40%) have predominantly disk morpholo- 
gies, inconsistent with major mergers fueling these sys- 
tems. In this case, the estimate of F from population 
statistics is not a true duty cycle, but is effectively di- 
luted by the population of non-merger QSOs. A third 
possibility is that unobscured activity may be substan- 
tially suppressed in the early stages of a merger. In such 
a scenario, merging SMBHs separated by more than a 
few kpc are much less likely to achieve high accretion 
rates and exhibit QSO luminosities than SMBHs close 
to or after final coalescence. This notion jibes with th e 
results of merger simulations (e.g.. iHopkins et "al1 l2008'). 
At present, given the limited statistics, we cannot easily 
distinguish which of these reasons is most significant in 
explaining the difference between the dual AGN fraction 
we measure and that expected from the statistics of the 
QSO population. All three may play a role to some de- 
gree. In addition, the dual AGN fractions we calculate 
here are applicable most directly only to moderate lu- 
minosity Type I AGNs. Dual fractions among luminous 
QSOs may be higher and systematic spectroscopic and 
imaging programs of luminous QSOs at high spatial res- 
olution could throw some light on the fueling properties 
of strongly accreting AGN. 

7. CONCLUSIONS AND FUTURE WORK 

From a well-defined sample of 12 SDSS radio-quiet 
Type I QSOs that display double-peaked structure in 
the [O III]A5007 line, we identify double structure in 6 
objects, using NIRC2 LGS AO imaging. We verify that 
both [O HI] velocity components in the QSOs are likely 
to come from AGN-ionized gas, consistent with the hy- 
pothesis that AGN exist in both galaxies in the six double 
systems. We model the images with GALFIT to estimate 
galaxy and point source luminosities, from which we esti- 
mate SMBH masses. There is no significant difference in 
Mbh between isolated and merging QSOs in our imaging 
sample, though the mean mass of our sample is higher 
than the parent SDSS population. Finally, we estimate 
that mergers detectable by current Keck LGS AO and 
NIRC2 imaging exist in 0.3%-0.65% of QSOs, which is 
considerably lower than the merger rate among inactive 
massive galaxies. 

In addition to developing improved morphology statis- 
tics via a larger imaging sample, an attractive next step 
would be to use AO spatially resolved IFU spectroscopy 
of our smallest-separation targets to answer the key ques- 
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tion of whether each galaxy in our pairs contains an ac- 
tive nucleus. We are currently pursuing this avenue are 
part of a follow-up study of the best dual AGN candi- 
dates. 

Data presented herein were obtained at the W. M. 
Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the 
University of California, and the National Aeronautics 
and Space Administration. The Observatory and the 
Keck II AO system were both made possible by the gen- 
erous financial support of the W. M. Keck Foundation. 
The authors wish to extend special thanks to those of 
Hawaiian ancestry, on whose sacred mountain we are 
privileged to be guests. Without their generous hospi- 
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